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Motivated by the recent measurement of the upper limit of B° —> J/ipD branching ratio, which is 
important in accounting for the soft J/ip production in B decays, we investigate B° —> J/ifD ^ and 
r/cD ^ decays in perturbative QCD approach based on kr factorization. Being pure annihilation 
(W-exchange) decays, these branching ratios are estimated to be at the order of 1CF 5 ~ 10“ 7 , 
which are just at the corner of being observable at the B factories. The measurements of these 
decay channels may help us to understand the QCD dynamics in the corresponding energy scale, 
especially the reliability of pQCD approach to these processes. 

PACS numbers: 13.25.Hw, 12.38.Bx 


I. INTRODUCRION 


In 1995, the CLEO Collaboration found a hump in the low momentum region of the inclusive spectrum 
of B —> J/if+ X decay Q. Later on, this observation was confirmed by Belle Q and BaBar Q. In these 
measurements, there is an excess in the momentum spectrum of the J/if recoiling mass at ~ 2 GeV. 
And, the excess corresponds to a branching ratio of 6 x 10 -4 . In order to explain this result, various 
hypotheses have been proposed 0,0,0|. 

In Ref. P|, Chang and Hou employ the idea of intrinsic charm jjj cc inside the B meson to this 
issue. Based on this scenario, they predicted that the branching ratio of B —> J/ifD should be about 
ICC 4 . However, according to recent BaBar and Belle measurements, the branching ratio upper limit of this 
process is less than 10 -5 jaBlj which implies that the intrinsic charm mechanism is not favored. In another 
scenario, in which the charmonium is produced predominantly in the Color-Octet mechanism, Eilam and 
Yang estimated the branching ratio of B —> J/ifD j|| and got a result of about 10 -8 . However, in the 
collinear factorization, they have to use a cut-off or (5-function to tame the end-point singularity. Hence, 
their numerical results are not stable. The recent progress in perturbative QCD (pQCD) treatment, 
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based on the fcy factorization, of B meson decays can solve this problem by introducing the Sudakov 


iy mi 

u 


form factor through the threshold resummation. Now, the pQCD approach |10| has become one of the 
broadly used theoretical methods in investigating the B meson two-body non-leptonic decays. Base on 


;ays. 

A 


the pQCD approach, many B meson decay modes have been calculated, like B —> Kir, mr llj, etc., and 
most results are consistent with the experimental data. Since there is no end-point singularity, the pQCD 


Q. 


approach can also be applied to the pure ’’annihilation processes ”, such as B —> D S K 

In this work, we calculate the B —> J/ipD M and B —> rj C D W processes in the pQCD kr factorization. 
In the decay of B —» J/ipD, the W boson exchange induces the four quark operator cb —> ud , and an 
additional pair of cc is created by a gluon. This gluon can attach to any quark involving in the four- 
quark operator. In the rest frame of B meson, the produced c and c quarks in the final states have the 
momenta of order 0{P^/ 2) and 0(Pd/ 2), respectively. Therefore, the gluon, which generates the charm 
quark pair, possesses a virtuality of order ~ 0{Mb/ 2), which enables the perturbative QCD calculation 
reliable. 

The paper is organized as follows: we present the formalism used in the calculation of B —♦ J/tpD^ 
and B —> r] c D ( P decays in Section El In Section EH we give out the numerical calculation results and 
some discussion on them. The last section is left for conclusions and summary. 


II. KINEMATICS 

The effective Hamiltonian for decay modes B —> J/ipD M and B —> rj c D^* ] is given by 



n cS 

Or 

o 2 



( 1 ) 

C7/x(! - 75)w d7 M (l - 75 )b, 


- 7 5 )u 07^(1 - 7 5 )b. 

( 2 ) 


As usual, in the pQCD approach the momenta of the final states are expressed in its light-cone compo¬ 
nents, like 


P = (P + ,P ,PT) = 


p u + pi pU _ pi 


V2 ’ V2 


,{p\p 2 ) 


( 3 ) 


And, the decay amplitude can be generally written as: 


M 


I dx\dx2dx3,bidbib2db2b2,db^ 


xTr[C(t)$B(xi,b 1 )$,p(x2,b2)$D(x3,b 3 )H(xi,b i ,t)e s(t) ] . 


( 4 ) 


Here, Tr denotes the trace over Dirac and color indices. C(t) is Wilson coefficient of the four quark 
operator which results from the radiative corrections at short distance. denote the wave functions 
which are process independent and represent the non-perturbative dynamics of hadronization. The hard 
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interaction kernel H is, nevertheless, process-dependent and can be calculated by perturbation QCD. t is 
chosen as the largest energy scale involving in the hard interaction to avoid the largest logarithms. S(t) 
is Sudakov form factor resulted from the resumnration of double logarithms [illfl3l . Therefore, in eq. 0 
only the hard part is process dependent and will be calculated in the following. 


A. The B —+ J/ipD Decays 


Of the B- and D W -meson wavefunctions, we make use of the same parameterizations as used in the 
s «ies of different _ Q Q Per vector meeon, in tortne of tire notation in Ref. Q. we 
decompose the nonlocal matrix elements for the longitudinally and transversely polarized J/ip mesons 
into 


(J/ip(P,eT)\c(z)jc(0)i\0) = 


1 


V2N~cJo 

1 f 1 


dxe lx ' z \m J/ T li y L \i j 'i> ( x) + [/ L P]i^ (x) ^ , 


dxe 


ixP-z 


mj/^yT\ij^ V (x) + [i-T P]ij^> T {x) 


( 5 ) 


( 6 ) 


V2N~cJo 

respectively. Here, and denote for the twist-2 distribution amplitudes, and T* and T 1 for the 
twist-3 distribution amplitudes, x represents the momentum fraction of the charm quark inside the 
charmonium. 

The J/'ip meson asymptotic distribution amplitudes read as [3] 


V L {x) = $ T (x) = 9.58 


2sfW c 


x{l — x) 


x(l — X ) 


^(x) = 10.94 ^-L (1-2x) 2 


T y (x) = 1.67 


2^2 N~ c 
fj/ ip 


l-2.8x{l-x) 
n 0.7 


0.7 


x(l — x) 


1 — 2.8x(l — x) 


[l + (2x — l) 2 ] 


x(l — x) 


-l 0.7 


1 - 2.8x(l - x) 


( 7 ) 


2^2 W c 

in which the twist-3 ones 'I> t v vanish, as the twist-2 ones, at the end points due to the factor [x(l — a:)] 0, '. 
In contrast to Ref. Q , here we distinguish the longitudinal and transverse distribution amplitudes of the 
polarized J/ip, which can exhibit the different asymptotic behaviors of these two types. 

From the effective Hamiltonian o, the Feynman diagrams corresponding to the concerned process are 
drawn in Fig.l, where the heavy dots denote the four quark operators. Similar figures can be obtained 
by replacing the J/ip by rj c for B —> rj c D process, and D by D* for the vector D -meson processes. With 
the meson wave functions and Sudakov factors, the hard amplitude for factorizable annihilation diagrams 
Figs. 1(a) and (b) is 

pi pOO 

F a = 16nC F Ml / dx 2 dx 3 / b 2 db 2 b 3 db 3 <p D (x 3 ) 

Jo Jo 

x (x 3 - 1 - x 3 rl - x 3 rl)^ L {x 2 )E f {t 1 )h 1 {x 2 ,x 3 ,b 2 ,b 3 ) - {[x 2 - 1+ 

(1 - 2 x 2 )rl + (1 - x 2 )r 2 ]'i’ L (x 2 ) + 2x 2 r 2 r 3 'b t (x 2 )}E f (t 2 )h 2 (x2,x 3 ,b 2 ,b 3 ) . (8) 
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FIG. 1: Feynman diagrams for B° —> J/4/D 0 decay process in pQCD. 


Here, the functions Ef(t a ) contain Sudakov factors and Wilson coefficients of four quark operator, and 
hard scale t a . The h a , the virtual quark and gluon propagator, are given in the appendix. 

The result for the non-factorizable annihilation processes, shown in Figs. 1(c) and (d), is 

M a = —JL=64 itCfMb / dx\dx 2 dx 3 / b\db\ b 2 db 2 </>b(xi, bi)(t) D {x 3 ) 

V J o J o 

x {(1 - 2 r\) (1 - x 3 ) ty L {x 2 ) +r 2 ( x 2 - x 3 ) r 3 ^> t (x 2 )}E m (t 3 )h 3 (xi, x 2 , x 3 , for, b 2 ) 

- {(1 - x 2 - (1 - 2x 2 )rl - (1 - 2x 2 + x 3 )rl)^ L (x 2 ) 

- r 2 (x 2 - x 3 )r 3 '& t (x 2 )}E rn (t 4 )h 4: (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) . (9) 

The total decay amplitude for this decay is: 


A a {B - J/1>D) = f B F a + M a . 


( 10 ) 


Thus, the B meson decay width of the concerned process is: 


T(B -> J/4>D) = 


g 2 f m% 

128t r 


(1 - r 2 - r 2 )\V cb V; d A a (B J/^D)Y 


( 11 ) 


B. The B — » rj c D^ Decays 


The nonlocal matrix element of r] c production from vacuum can be generally expressed as 

{Vc(P)\c(z)jc{0)i\0) = ~^= dxe wP - z ^[) 5 P]ijrf{x) + j . (12) 
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Here, r) v (x ) and r/ s (x) denote the twist-2 and twist-3 r] c meson distribution amplitudes, respectively. The 
asymptotic forms of the ?y c distribution amplitudes are given in |l?| | : 


r] v {x) = 9.58 


77 s (a;) = 1.97 


fv 


2 ^ 2 TVc 

fvc 


x{\ — x) 


x(l — x) 


-t 0.7 


1 - 2.8z(l - x) 


7(1 — x) 


0.7 


,_ , x . (13) 

2y/2Nc [l-2.8x(l-a;)_ v ' 

Performing the similar procedure as in above subsection, we can get the decay amplitudes for B —> r] c D 

and B —> r) c D* straightforwardly. 


C. The B —> J/ipD* Decays 


The B —> J/ipD* decay rate are 


r = 


g 2 f p c 

32t tMI 

B a—L,T 


, 


(14) 


where P c = \P 2z \ = |P 3 Z | are the momenta of the outgoing vector mesons; the superscript er denotes 
for the helicity states of the two vector mesons, the L for the longitudinal and T for the transverse 
components. The amplitude can be decomposed, according to the Lorentz structure, to fist : 


-4 (CT) = 


«r + 


b c 

-p»pv + i. 


A va0 P 2a Pi8 


m^rriD m^mn 

= MlA L + M 2 B A N e* 2 (a = T ) • e*(a = T) + iA T e a ^ p e* 2a (a)e* 30 (a)P 2l P 3p , 

with the convention e 0123 = 1 for the total anti-symmetric tensor and definitions 


M 2 B A l = a e 2 (L) ■ e 3 (L) + 

Ml A N = a e* 2 (T) ■ e*(T) 
c 


m^rriD 


e* 2 (L)-P 1 et(L)-P 1 , 


At 


1 ;■ 

Hereby, the only work left is to calculate the matrix elements Al, An and At with 

A i = f B F i +M i ,(i = L,N,T). 

Here, i 7 ) and Mi, coming from the calculation of hard interaction, are given as follows: 

pi pOO 

F L = IQiyGfM'I / dx 2 dx 3 / b 2 db 2 b 3 db 3 4>d(x 3 ) 

Jo Jo 

x {x 3 1 (x 3 - 2)7-3 - ( 2x 3 - l)rl)^ L {x2)Ef(t 1 )h 1 {x2,x 3 ,b 2 ,b 3 ) 

+ {{x 2 - 1 + (1 - 2 x 2 )r 3 + (1 - X2)rl]^ L (x2)}E f (t2)h 2 (x2,x 3 ,b2,b 3 ) 


(15) 


(16) 


(17) 


(18) 


pi pOO 

F n = IMC F M 2 B / dx 2 dx 3 r 2 r 3 / b 2 db 2 b 3 db 3 <f)D{x 3 )r 2 r 3 
Jo Jo 

X {l-x 3 )'S/ v (x 2 )Ef(t 1 )h 1 (x 2 ,x 3 ,b 2 ,b 3 ) + {(x 2 -2)'l! v (x2)}Ef(t 2 )h2(x 2 ,x 3 ,b 2 ,b 3 ) , (19) 
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pi /»00 

F T = 32nC F Mg / dx 2 dx 3 r 2 r 3 / b 2 db 2 b 3 db 3 <j) D (x 3 )r 2 r 3 
Jo Jo 


(1 + X 3 )'I' V ( x 2 )Ef (ti)hi(x 2 , x 3 , b 2l b 3 ) + x 2 ty v (x 2 )Ef(t 2 )h 2 (x 2 ,x 3 , b 2 , b 3 ) , (20) 

M l = —J=MtiC f M% f dx\dx 2 dx 3 f b\db\ b 2 db 2 <f> B (xi, b\)(j> D {x 3 ) 

\2JM C Jo Jo 

x {(cc 3 — 1 — 2(x 3 - l)(r| + r 2 )) T L (x 2 ) +r 2 (x 2 + x 3 - 2)r 3 ^> t (x 2 )} 
x E m (t 3 )h 3 (xi,x 2 ,x 3 ,b 1 ,b 2 ) + {(1 - x 2 - (1 - 2 x 2 )r% - (1 - 2x 2 - x 3 )rl)'S> L (x 2 ) 

- r 2 (x 2 + x 3 )r 3 ^ t (x 2 )}E rn (t 4 )h 4: {xi,x 2 ,x 3 ,bi,b 2 ) , (21) 

1 f 1 

Mjy = 647 tC f M f / dx\dx 2 dx 3 / b\db\ b 2 db 2 <f) B (x\, b\)<j)D{x 3 ) 

v2A c Jo Jo 

x {(^3 - 1) r 3 'S> T (x 2 ) + (x 2 - 1) r%'P T (x 2 )}E m (t 3 )h 3 (x 1 ,x 2 , x 3 , b 1 ,b 2 ) 

+ {(-x 3 rl^ T (x 2 ) + 2r 2 r 3 ’$’ v (x 2 ) - x 2 rj'i> T (x 2 )}E m (t 4: )h 4 (x 1 ,x 2 ,x 3 ,bi,b 2 ) , (22) 

1 9 f 1 f°° 

M T = 12 HttC f M^. / dx\dx 2 dx 3 / bidbi b- 2 db 2 <f) B (xi, b\)4>D{x 3 ) 

v2 N c Jo Jo 

x {(x 3 - l)r^ T (x 2 ) - (x 2 - l)r 2 'i! T (x 2 )}E m (t 3 )h 3 (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) 

+ {(-x 3 rl'& T (x 2 ) + x 2 rl'5> T (x 2 )}E m (t4)h4(xi,x 2 ,x 3 ,bi,b2) ■ (23) 


III. NUMERICAL RESULTS 


In this work, the input parameters for the numerical calculation are 
literature, 



, which are commonly used in 


m JN , = 3.097 GeV, m Vc = 2.980 GeV,/ D . = 230 MeV, 
f D = 240 MeV, = 1.87 GeV, m D * = 2.005 GeV, 
m B = 5.28 GeV, \V cb \ = 0.043, \V ud \ = 0.975, t b = 1.54 x 10” 12 s. (24) 

At leading order, the main uncertainty comes from the meson wave functions. Fortunately, the meson 
wave function, that describes hadronic process, is universal at a certain scale. For instance, the B meson 
wave function is constrained by the measured exclusive hadronic decays, like B — > 7T7t. tt IhI with 
parameter uj b from 0.32 to 0.48. To determine the D meson wave function is more tough task than that 
of B meson, because the heavy quark limit here is not as good as in the B meson case. Referring to to 
B -> DMmQ process, we can fit the D meson wave function parameter to be a B = 0.8 ± 0.2. The 
charmonium distribution amplitudes can be inferred from the non-relativistic heavy quarkonium bound 
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state wave functions, which have been shown to be successful in describing the charmonium production 


in e + e _ collisionsThe meson decay constant can be measured via its pure leptonic decay. We have 
fj/,p = 405 ± 14 MeV and f Vc = 420 ± 50 MeV. In addition to the uncertainties remaining in the above 
input parameters, the higher order corrections to the hard part are also important, which is discussed in 
Ref. l2o|. 

Considering of the above uncertainties discussed, we can give out the branching ratios of the discussed 
processes with error bars: 


Br(R° -> J/ipD) 
Br (B° -> Vc D) 
Br (B° -> tj c D*) 
Br (B° -> J/ipD*) 


(3.45t};^ ± 1-51 ± 0.32) x 1 (T 6 , 
(1.28l^i ± 0.58 ± 0.35) x 1(T 5 , 
( 8 . 26 ^ 2.34 ± 2.23 ± 2.06) x 1(T 6 , 
(7.04^2 54 ± 2.72 ± 0.53) x 1CT 7 . 


(25) 


In the above, the uncertainties mainly come from u>b, o-d, and the decay constants, respectively. To 
diminish the uncertainties, for B° —> J/ipD* process, we evaluate the longitudinal polarization fraction, 
that is: 

Pl = ^ = 0 . 66 . (26) 

This polarization fraction is not sensitive to the above mentioned input parameters, because they only 
give an equally change of each polarization amplitudes. However, this fraction is still sensitive to the J/ip 
wave function. If we set the distribution amplitude of transversal part the same as longitudinal part, the 
branching ratio become larger and the polarization fractions changed: 

Br(B° -> J/i/jD*) = 10.5 x 10~ 7 . (27) 

P L = ^ = 0.40; P n ,t = ^ = 0.30. (28) 


That is to say that for B° —> J/ipD* the most important uncertainty comes from the vector meson wave 
functions. 


f. B, 


Compared to Ref. 6 |, our results are much bigger. In | 6 |, all wave functions, which describe the non- 
perturbative hadronization, are 5-function-like. However, the 5-like wave function can not embody the 
relativistic corrections, though it can be used to avoid the end-point singularity due to the wave function 
overlap absent. In this work, the hadron distribution amplitudes are obtained from from the established 
models with experimental fittings. In our work, we take into account the Sudakov form factor and the 
transverse momentum kj* distribution, which are unique characters of pQCD approach. For B° —> J/ipD* 
process, since the charmonium longitudinal distribution amplitude is different from its transverse one, 
and hence our longitudinal polarization fraction are larger than what obtained in Ref. jfj]. 

Since there is only one kind of CKM phase involving in the concerned process, there should be no CP 
violation in these process within the standard model. On experimental side, so far there is only an upper 



limit for the branching ratio of B° —> J/%l)D process. That is 


Br {B° -► J/i>D) < 1.3 x 10“ 5 [8], 

Br (B° -► J/i>D) < 2.0 x 1CT 5 [9], 

from different experiment group, which is larger than, but very close to our prediction. 


( 29 ) 


IV. SUMMARY 

In this work, we have calculated the decays of B° —> J/tp(j] c )D^ in the pQCD approach. These B 
meson exclusive decay processes are in pure annihilation type, which is hard to be accurately calculated 
in other approaches with the end-point singularity. By keeping the transverse momentum hr, the end¬ 
point singularity disappears in our calculation. Our numerical results shows that the branching ratios 
of B° —> r) c D, B° —> r] c D* , B° —> J/i)>D and B° —» J/ipD^ decay processes are of the order 10 -5 , 
10 -6 , 10~ 6 , and 10 -7 , respectively, which is just close to the experiment capability to measure them. 
Although both Belle and BaBar measured the J/ip momentum spectrum in B inclusive decays, they 
did not obtain the branching ratios of these exclusive decays modes. Considering that the upper limits 
set by experiments are very close to our predictions. We suggest that BaBar and Belle measure these 
exclusive processes in near future. The observation of these exclusive processes may greatly improve our 
understanding on the B meson exclusive hadronic decays, and the corresponding theory describing them 
as well. 
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APPENDIX A: SOME FUNCTIONS 

The function E l f, E m , and E' m including Wilson coefficients are defined as 

E f (t) = (Ci(t) + ^j^)a s (t) 


(Al) 


E m (t ) = C 2 (t)a a (t) e 


— S B{t) — S<b (t) — S D{t) 


(A2) 
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where S B , SW, and Sd result from summing both double logarithms caused by soft gluon corrections and 
single ones due to the renormalization of ultra-violet divergence. The above Sb,^,d are defined as 


S B (t) = s(xiPf ,bi) + 2 f 

Jl/bi h 

S*(t) = s(x 2 P^,b 3 ) + 2 f ^77g(M , )> 
J i / 6 2 F 

S D {t) = s{x 3 P 3 ,b 3 ) + 2 f ^-7g(/x'). 

Ji/b 3 A* 


(A3) 

(A4) 

(A5) 


where s(Q,b), so-called Sudakov factor, is given in Reference [ 21 ] . 

The functions hi— 1 , 2 , 3,4 hr the decay amplitudes come from the propagator of virtual quark and gluon. 
They are defined by 

hi(x 2 ,x 3 ,b 2 ,b 3 ) = (x 2 - l)(x 2 - x 3 )rl - (x 3 - l){x 2 r 3 - x 2 + l)b 2 ) 

x - x 3 + x 3 r% - r| b 2 )J 0 (M B \J 1 - x 3 + x 3 r% - rf b 3 ) 6 (b 2 - b 3 ) + (b 2 <-> 63)|, (A6) 

h 2 (x 2 ,x 3 ,b 2 ,b 3 ) = H^\m B \J ( x 2 - l)(x 2 - x 3 )r% - ( x 3 - l)(a; 2 r§ - x 2 + 1) b 2 ) 


x - x 2 + x 2 r\ - x 2 r\ b 2 )J 0 (M B yJl - x 2 + x 2 r| - x 2 rj b 3 )0(b 2 - b 3 ) + (b 2 ^ & 3 ) j 


(A7) 


h(j= 3t 4 )(xi,x 2 ,x 3 ,bi,b 2 ) = 

ttH (x 2 - l)(x 2 - x 3 )r\ - (x 3 - l)(x 2 r\ -x 2 + l) h) 

x J 0 (M B \J (x 2 - l)(x 2 - x 3 )rl - (x 3 - 1 )(x 2 r§ - x 2 + 1) b 2 )9(bi - b 2 ) 


+ (bi ^ b 2 )\ x 


K 0 (M B F {j) b r), for Ffa > 0 

^\Mb J\F^\ h), for Ffa < Oy 


(AS) 


where Hq 1 ^) = Jq(z) + i Y 0 (z), and Fij \s are defined by 


^( 3 ) = {xx+x 2 -l)(x 2 -x 3 )rl + (x 3 -l)(xi+x 2 -l-x 2 rl), 
F ( 4 ) = -(aq - x 2 )(x 2 - x 3 )r% - x 2 x 3 r\ - x^x 3 + x 2 x 3 - 1. 


(A9) 


The hard scale t’s in the amplitudes are taken as the largest energy scale in the FI to kill the large 
logarithmic radiative corrections: 


t\ = m.ax(M B \J 1 - x 3 + x 3 r% - r\ , 1 /b 2 , l/b 3 ), (A10) 

t 2 = ma x(M b J 1 - x 2 + x 2 r% - x 2 r|, l/b 2 , l/& 3 ), (All) 

i? = max.(M B ^l\F^\,M B \J (x 2 - l)(x 2 - £ 3 )r| - (a; 3 - 1 )(x 2 r§ - x 2 + 1), 1/for, l/fo 2 )- (A12) 
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